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Bulk QCD thermodynamics and sterile neutrino dark matter 
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We point out that the relic densities of singlet (sterile) neutrinos of interest in viable warm and 
cold dark matter scenarios, depend on the characteristics of the QCD transition in the early uni- 
verse. In the most promising of these dark matter scenarios the production of the singlets occurs 
at or near the QCD transition. Since production of the singlets, their dilution, and the disappear- 
ance of weak scatterers occur simultaneously, we calculate these processes contemporaneously to 
obtain accurate predictions of relic sterile neutrino mass density. Therefore, a determination of 
the mass and superweak mixing of the singlet neutrino through, for example, its radiative decay, 
along with a determination of its contribution to the critical density, can provide insight into the 
finite-temperature QCD transition. 
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I. INTRODUCTION 

In this paper we show how the temperature and the 
order (first order or not) of the QCD transition in the 
early universe is important in determining the relic densi- 
ties of sterile neutrinos in models where they provide the 
dark matter. The QCD transition epoch is the regime 
where the quark and gluon degrees of freedom annihi- 
late, chiral symmetry is broken, and quarks are incorpo- 
rated into color singlets. For the low chemical potential 
(i.e., small baryon number or high entropy) characteriz- 
ing the universe, we expect the QCD epoch to occur at 
a temperature (energy scale) of order of the pion mass, 
T cr it ~ 100 MeV. Coincidcntally, this is the tempera- 
ture regime where the dark matter singlet neutrinos are 
produced by scattering-related processes. 

Therefore, an accurate calculation of the density of 
dark matter candidate sterile neutrinos requires detailed 
specification of the evolution of the background plasma 
through the QCD transition. It is conceivable that this 
relation between the QCD transition and sterile neutrino 
properties could be exploited to understand the cosmic 
QCD epoch. This approach is not available, for exam- 
ple, with supersymmetric dark matter candidates which 
freeze-out of the plasma at temperatures well above the 
QCD transition. Here, we also review the current best 
constraints on sterile neutrino dark matter candidates. 

The current cosmological paradigm of a spatially flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) uni- 
verse with critical density components shared between 
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cold dark matter [(CDM), f2 m ~ 0.3] and a negative- 
pressure dark energy (Q x ~ 0.7), with nearly scale- 
invariant adiabatic Gaussian perturbations produced by 
inflation leading to structure formation (ACDM), is sup- 
ported by a concordance of observational evidence from 
the cosmic microwave background, large scale structure, 
rich galaxy clusters, and high-redshift Type-la super- 
novae |Q. There is much optimism in observational cos- 
mology, but a number of unresolved problems remain. 
Among them is the very nature of particle dark matter 
and dark energy. 

More astrophysically, there may be a discrepancy be- 
tween the quantity of substructures on galactic scales pre- 
dicted in numerical simulations using ACDM initial con- 
ditions and that observed in the local group: the "dwarf- 
galaxy problem" ^ [|. Another discrepancy exists be- 
tween the observed flat and low density profiles of dwarf 
galaxy cores M || relative to those found in A-body 
simulations. One solution to the dwarf-galaxy problem 
is a modification of the primordial power spectrum at 
sub-galaxy scales by thermal damping of perturbations 
by warm dark matter (WDM) [0. Dwarf-galaxy halos 
in WDM scenarios form through fragmentation of larger 
structures, and the concentrations of their halos may be 
reduced 0. Although these problems may be signaling 
the presence of WDM, the apparent dwarf-galaxy prob- 
lem may be the result of dispersion of baryons from halos 
due to late reionization or supernovae, so that the unob- 
served dwarf halos are simply dark. These issues could be 
resolved through refined calculations of star formation in 
small halos, or through observations of galactic substruc- 
ture in, for example, gravitational lensing studies [||. 

Sterile neutrinos are a natural candidate for WDM || 
|l0| and CDM (TTJ, [T^j , a nd can emerge from models with 
composite fermions floj] , mirror fermions |l4j or light ax- 
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inos |LE|. Sterile neutrinos with masses of 1 — 100 keV 
and small mixings with active neutrinos are produced 
in a nonequilibrium quantum decoherence process as- 
sociated with the scattering of active neutrinos: either 
with the usual simplifying assumption of a small univer- 
sal lepton number || (of order of the baryon number) or 
with the more liberal possibility of a large lepton num- 
ber in neutrinos (several orders of magnitude larger than 
the baryon number) pj] |. We define the lepton num- 
ber in a neutrino flavor here as the difference between 
neutrino and antineutrino number densities normalized 
by the photon number density, n 7 , at a given epoch: 
L a = (n Ua - np a )/n*y. 

Constraints on the lepton numbers which reside in the 
neutrino seas can be made by requiring successful pro- 
duction of the light element abundances in primordial 
nucleosynthesis while keeping the baryon density com- 
patible with that inferred from observations of the cos- 
mic microwave anisotropics Jl6| . Because of the near 
bi-maximal mixing paradigm emerging for the neutrino 
mass matrix, more stringent constraints can be placed 
on the neutrino lepton number since asymmetries in the 
more poorly constrained mu and tau flavors can be con- 
verted to an asymmetry in the more stringently con- 
strained electron neutrinos through synchronized flavor 
transformation [i7j |. However, even the strongest con- 
straint available, on the electron neutrino lepton num- 
ber, allows |i e |iMeV 0.1 at the epoch of nucleosynthe- 
sis. Because sterile neutrinos are produced at temper- 
atures much higher than those characterizing big-bang 
nucleosynthesis (BBN), and because the lepton number 
(unlike degeneracy parameter) is not a comoving invari- 
ant, we quote values of the lepton number at the epoch 
where T = 100 GeV. Therefore, the most conservative 
upper bound on the lepton number at this temperature is 
weaker by the ratio of the number of statistical degrees of 
freedom g(T) in the plasma, #(100 GeV)/g(l MeV), or 
I -f^ | ioo GeV ^5 1- All lepton- number-driven mechanisms 
we consider are well within this conservative limit. 

Depending on their masses and non-thermal energy 
distributions, the sterile neutrinos produced in the early 
universe can be CDM, WDM, or hot dark matter (HDM). 
The production mechanisms and constraints on these 
models were studied in Rcfs. [TL0, 12 . The strongest con- 



straints on the upper bound of the mass of the sterile 
neutrino derive from x-ray emission from nearby galaxy 
clusters, viz. the Virgo cluster m s < 5 keV when 
L ~ 10~ 10 . Lower bounds on the mass of the WDM par- 
ticle can be made by requiring that its free-streaming not 
excessively damp out the small-scale structure, i.e., be- 
have too much like HDM. Limits of this type come from 
the observed structure on small scales in the high-redshift 
Lyman-a forest ||l9|| and by requiring sufficiently early 
rcionization of the universe seen by the Gunn-Peterson 
effect [^0|. A lower bound on the mass can also be de- 
rived from the observed phase-space densities of galaxy 
cores [§). All of these limits require out-of-equilibrium 
("standard") WDM particles with mass m > 0.75 keV, 



or sterile neutrino WDM in models where L ~ 10~ 10 
with mass m s > 2.6 keV ^T|. It is exciting that the cur- 
rently unconstrained mass range 2.6 keV < m s < 5 keV 
(at L ~ 10~ 10 ) is that favored by WDM solutions to 
small scale structure formation, and, even more, this 
mass range for the sterile neutrino dark matter is also 
eminently detectable via radiative decay p8[ . 



II. THE QUARK-HADRON TRANSITION 

The study of the condensed matter physics of QCD 
remains an active topic, through analysis of the full color 
SU(3) QCD Lagrangian with numerical lattice simula- 
tion and analogy with spin systems (for a recent survey 
see Ref. p2[). The nature of the QCD transition at finite 
temperature and zero chemical potential /i q is clear in 
models containing two massless quarks (m qu — m qd = 0, 
m qB —>■ oo): the transition is second order With 



physically realistic small but nonzero m 9u and 



the 

second order transition loses its critical behavior and is 
a smooth crossover |24|] . In degenerate three-flavor mod- 
els (m qu = m qd = m qs = 0), it was inferred analyti- 
cally that the transition is first order and this has 
been supported by lattice studies p5[ . Therefore, with 
fi q — * (as is appropriate for the early universe where 
(iq/T ~ 10~ 8 ) there exists some critical mass m qa , below 
which the transition is first order. The position of this 
critical point in temperature, and m qg is still under 
investigation by lattice simulations and may be probed 
by relativistic heavy ion collision experiments p4j . Lat- 
tice calculations indicate m qs is well below the inferred 
physical strange quark mass |25|, ^6|, but these results 
are not conclusive p7| . Since production of sterile neutri- 
nos peaks at T ~ 130 MeV(m s /3 keV) 1 / 3 , and the QCD 
transition likely occurs from T cr i t = 100 — 200 MeV, the 
effects of the transition on production are unmistakably 
important. 

Convincing theoretical motivation or experimental de- 
tection of the critical end point as yet does not exist, 
and so whether the QCD transition in the early universe 
behaves as condensing water (first order) or thickening 
pudding (smooth crossover) remains to be discovered. 
We argue below how the dark matter, if sterile neutri- 
nos, can also be used to probe the critical behavior of 
QCD. 

We model the evolution of the universe through a first- 
order transition as in Refs. |^, The transition from 
the high-temperature quark-gluon phase to the hadron 
phase includes the change in the statistical degrees of 
freedom of the plasma, the release of vacuum energy, and 
the annihilation of quark- antiquark pairs. Supercooling 
of the plasma can lead to nucleation of remaining quarks 
into dense isothermal baryon number fluctuations [j30| , 
but this does not affect the analysis here. 

A first-order transition evolves through a constant- 
temperature epoch, with a duration of the order of a 
Hubble time at that epoch. During this constant tern- 
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perature, mixed phase regime, volume is swapped from 
the quark-gluon phase to bubbles of the hadronic phase 
as the universe expands. In this scenario the fraction of 
the volume in the quark-gluon phase evolves from f v — 1 
to 0, at which time the transition ends. The two phases 
differ in entropy density and therefore latent heat is re- 
leased in the transition. The dynamics of the scale factor, 
a, through the transition is given by (2^, |3l|| 
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where x = g q /gh is the ratio of the number of relativistic 
degrees of freedom in the quark g q and hadron g^ phases 
(approximately valid during the mixed phase epoch), and 
the characteristic QCD expansion rate is 
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Here the vacuum energy is parametrized by an effective 
"bag constant" B. Pressure equilibrium between quark 
and hadron phases is used to find the relation between the 
vacuum energy and the transition's critical temperature, 

-^crit ■ 

Using the conservation of comoving entropy it can be 
shown that the fraction of the plasma in the quark-gluon 
phase at a given time, t, from the start of the transition 
at time, U, is EQl 
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These relations are used to calculate the cosmological dy- 
namics for sterile neutrino dark matter production dur- 
ing the first-order transition. The evolution of tempera- 
ture for a first-order and crossover transition is shown in 
Fig. | 



III. DARK MATTER PRODUCTION 

Production of sterile neutrino dark matter occurs 
through decoherence of the active neutrino gas via colli- 
sions with the background plasma of weakly interacting 
particles. We calculate this nonequilibrium production 
of sterile neutrinos as described in Rcf. [Q, which we 
summarize here. This production is described by a semi- 
classical Boltzmann equation that captures the quantum 
behavior of the system. The distribution function for 
sterile neutrinos f s over momenta p evolves as 

|/.(p,t)-ffp|/.(p,*) 



T(v a ->■ v s ;p,t) [f a (p,t) - f s (p,t)} 



(3.1) 
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FIG. 1: Evolution of temperature for a first-order transition 
(solid) and crossover transition (dashed) for T cr i t = 150 MeV. 



Here f a is the active neutrino distribution function that 
is coupled to the sterile, T(v a — > v s \p,t) is the effective 
rate of production (annihilation) of sterile states and H = 
a/a. 

It is important to note that the scattering kernels 
T(y a — > v s ',P, t) for the case of sterile neutrino production 
are always momentum conserving [j$2) so that Eq. (3.1) 
is exact. We handle the calculation of this scattering ker- 
nel as described in Ref. [ fl2[ . This includes augmenting 
the scattering rate as the population of weak scatterers 
increases at high temperatures. 

The behavior of the QCD transition affects dark mat- 
ter production through the varied evolution of the total 
energy density, ptot, of the plasma for different critical 
behavior. In turn, this modifies the evolution of the 
expansion rate through the Friedmann equation H = 



y/&K/3r 



i 1/2 

rn pl p tot . The dynamics of expansion enters the 
calculation of dark matter production in two ways: (1) 
alteration of the time-temperature relationship; (2) the 
modification of redshifting of relativistic species, includ- 
ing sterile neutrinos. 

If the universe underwent a crossover QCD transition 
(i.e., no phase transition), the effects on sterile neutrino 
densities stem entirely from the annihilation of quark- 
antiquark pairs and the resulting disappearance of de- 
grees of freedom in the plasma. The entropy carried in 
the annihilating species is transferred to the remaining 
constituents of the plasma. The time-temperature rela- 
tion in the early universe with evolving numbers of ther- 
malized species is far from a new problem. This is basi- 
cally the same effect (though more pronounced) on the 
time-temperature relation as that caused by the annihi- 
lation of electron-positron pairs when T ~ 25 keV. A 
general treatment of the time-temperature relationship 
including the effects of the population of sterile neutrinos, 
is given in the Appendix of Ref. |l2| , whose approach we 
have incorporated. A crossover transition disturbs pro- 
duction by causing the universe to "linger" for some time 
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FIG. 2: The fraction of critical density Q a produced above 
a temperature r pro( j. Here we have taken T CT i t = 150 MeV. 
There is significant enhancement due to thermal lingering af- 
ter T cr it for a crossover transition (dashed) , and there is a dra- 
matic enhancement at T cr it in the first order case (solid) . The 
most significant difference between the two cases stems from 
the delay of the radiation-dominated era in a first-order tran- 
sition. Here we have taken m s — 3 keV, sin 2 26 = 2.6 x 1CP 8 , 
and L ~ lO -10 . 



at temperatures just below the critical temperature (see 
Fig. [|) . Since the QCD transition temperatures can be at 
or near the peak rate of dark matter production, the po- 
sition of the post-critical temperature lingering may en- 
hance the production. A representative case of dark mat- 
ter production through a crossover transition is shown in 

Fig. |. 

The temperature is constant in a first-order transition, 
and the duration of the transition is given by the time re- 
quired for the quark phase to disappear, Eq. (2.2). Here, 
the production with temperature is a step-like function 
at the critical temperature as a result of the finite time 
spent at T cr ;t (see Fig. ^). More significant than the pro- 
duction during the transition is the delay of the radiation 
domination epoch by a first-order transition, so that the 
time-temperature relation is modified from the standard 
case, t oc T~ 2 . The duration of the first-order transition 
can be obtained by Eq. (O), 
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ferent transition scenarios. The expansion also modifies 
the redshifting of sterile neutrino momenta (poca -1 ). In 
a radiation-dominated environment, and with a crossover 
transition, p oc i" 1 / 2 . However, during a first-order phase 
transition the momentum will scale inv erse ly with the 
scale factor as found by integrating Eq. (2.1). Handling 
the redshifting explicitly through these scalings as op- 
posed to using the second term on the left-hand side of 
Eq. (3.1) greatly accelerates numerical computation by 
eliminating the calculation of the momentum derivative 
of the time-evolving sterile neutrino distribution. In this 
way, the different dependence of the scale factor with 
time gives rise to different redshift histories in different 
transition scenarios. 

The final effect of the transition we consider is sterile 
neutrino dilution by the heating of the thermally coupled 
species in the plasma relative to the uncoupled dark mat- 
ter by entropy transfer (latent heat release) . Exclusive of 
all other effects in a crossover transition, this produces a 
heating of the coupled species relative to the dark mat- 
ter. With the assumption that the crossover transition 
is rapid relative to the expansion time, the dilution of 
the dark matter relative to coupled species is given by 
comoving entropy conservation as x^ 1 = ghjg q ~ 1/3. 
A crossover transition does not necessarily take place in 
less than the expansion time scale, as can be seen for the 
cases we consider in Fig. |[ Therefore, the dilution should 
be handled by a Boltzmann formula tion as is automat- 
ically done in our treatment of Eq. (3T) and the time- 
temperature relation. For a crossover transition, dilution 
is the major factor in determining the different critical 
density contributions shown in Fig. |3|. In a first-order 
transition, dilution effects are handled in the evolution 



of the scale factor (2.1) and time-temperature relation. 
We have included all of the above effects in a numerical 



solution of the production equations (3.1). Contours of 
predicted sterile neutrino critical densities of Vl s h 2 = 0.15 
for three representative lepton number cosmologies and 
both first-order and crossover QCD transitions are shown 
in Fig. U The shape of the large lepton number cos- 
mology contours can be understood from the resonance 
condition for a neutrino of momentum p: 
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which for typical parameters is approximately 2 x 10 s. 
A first-order transition universe will then spend slightly 
more time reaching a given temperature. This can be 
seen in Fig. |l|. This delay in the temperature decrease 
with expansion is a major factor in enhancing dark mat- 
ter production in a first-order transition universe. That 
is, the partial equilibration of the sterile neutrinos has 
more time to operate in the first-order case. 

In this way, the different dependence of the scale factor 
with time gives rise to different redshift histories in dif- 



The largest disparity between first-order and crossover 
transitions occurs when the resonance position is in a 
populated portion of the neutrino momentum distribu- 
tion (of order p/T ~ 1) at T cr i t . For an L — 0.01 cos- 
mology, the resonance is in the peak of the momentum 
distribution at T cr j t for masses m s rs 13 keV, the pro- 
nounced effects of the nature of the transition for this 
case can be seen in Fig. [I| Furthermore, for this lepton 
number, the resonance position falls at very low momenta 
for m s < 1 keV during the peak production. Therefore, 
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FIG. 3: Fraction of critical density contributed by relic ster- 
ile neutrinos for a varying critical temperature for the QCD 
transition. Solid line is for a simple crossover QCD transi- 
tion; dashed line is for a first-order transition. Here we have 
specified m s = 3 keV, sin 2 29 = 2.6 x 1(T 8 , and L ~ 1(T 10 . 
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FIG. 4: Shown is the parameter space available for sterile 
neutrino dark matter, with varying lepton-number cosmolo- 
gies. The contours (numbered by their initial lepton number) 
are positions in the mass and mixing angle space where ster- 
ile neutrinos produce critical densities of Q B h 2 — 0.15. The 
thin (thick) lines are for first-order (crossover) QCD transi- 
tions (T cr it = 150 MeV). Also shown are the excluded re- 
gions (shaded gray) from small scale structures — the Gunn- 
Peterson bound and Lyman-a forest — and halo phase space 
densities, the resolution of the diffuse x-ray background by 
Chandra, and observations of the Virgo cluster by XMM- 
Newton. The dashed region is that which may be probed 
by the proposed Constellation- X mission ]18|. 



production is thermally suppressed, leading to the break 
seen at m s ~ 1 keV for L = 0.01 cosmologies in Fig. [|. 

An effect present above the quark-hadron transition is 
that of quark screening, or many-body neutrino-quark 
scattering. The mean distance between quarks above 
the quark-hadron transition at T ~ 150 MeV is approx- 
imately I ~ [3C(3)g,/7r 2 ]- 1 / 3 T- 1 - 1 fm, and the de 



Broglie wavelength of a typical neutrino is A ~ Ti/p ~ 
0.4 fm. Therefore, neutrinos below the average momen- 
tum will have A > I, and two body single-particle scat- 
tering is no longer a complete description. Many-body 
coherent scattering will tend to cancel the weak charge 
of quarks and antiquarks, and reduce the effectiveness 
of decohering interactions of active neutrinos. The ef- 
fects of many-body scattering are beyond the scope of 
this work. Though it can be significant in altering scat- 
tering rates above the QCD transition, the majority of 
sterile neutrino dark matter is produced below the QCD 
transition (as seen in Fig. ||) and the resulting effects of 
quark screening will be relatively small. 



IV. PROBING THE QCD TRANSITION WITH 
DARK MATTER 



In order to predict the sterile neutrino dark matter 
density one must specify the singlet neutrino rest mass 
m s , vacuum mixing with active neutrinos sin 2 29, the ini- 
tial lepton number L, and the order of the QCD transi- 
tion and its critical temperature T cr j t . Since the dark 
matter density is being determined rather precisely [j33| , 
the sterile neutrino production mechanism will constrain 
the relationship between the sterile neutrino properties 
and parameters describing the background plasma. 

It was shown in Ref. Jig] that the radiative decay of 
sterile neutrino dark matter v s — > v a ^ may be detected 
with high-sensitivity spectrographs aboard the current 
Chandra or XMM-Newton x-ray telescopes, or with the 
proposed large-surface-area and precision spectrograph 
of the Constellation X mission. As noted earlier, the 
lack of a significant line in the observation of the Virgo 
cluster constrains the largest possible mass of a sterile 
neutrino dark matter candidate produced in the zero- 
lepton- number case (m s ^ 5 keV) . This constraint comes 
from the fact that the Virgo cluster is one of the better 
candidate objects fl8|| and a high-sensitivity spectrum is 
available |34j . Observations of this and other clusters and 
field galaxies with current x-ray observatories can cither 
detect or exclude sterile neutrino dark matter as a dark 
matter candidate in zero-lepton-number universes. 

If sterile neutrinos are the dark matter, detection of ra- 
diative decay by a sufficiently sensitive observatory would 
readily identify m s since the radiative decay produces 
monoenergetic photons with energy E 7 — m s /2. Detec- 
tion of these photons in the x-ray is most likely in long- 
duration exposures of nearby high surface dark matter 
mass density objects. The x-ray flux coming from the 
radiative decay of a sterile neutrino is 
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Note that the radiative decay rate of a heavy sterile neu 
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trino is greatly enhanced due to the lack of the Glashow- 
Iliopoulos-Maiani suppression present in the decay of ac- 
tive neutrinos of the same mass and mixing ||^] . 

Since the position of the x-ray line will specify m s , the 
strength of the line determines sin 2 29, and therefore the 
position in Fig. |J. Of the remaining parameters (L a , 
T cr it and the order of the transition), only a small range 
will successfully produce the observed critical density and 
detected radiative decay flux. However, there remains a 
degeneracy between L a and the order of the transition 
since both a first-order transition and a slightly larger 
L a can enhance production at smaller mixing angles (see 
Fig. ||). However, if one can independently determine 
either the order of the transition or L a , then the observed 



flux [Eq. (4.1)1 w ill determine the remaining degree of 
freedom. 



V. CONCLUSIONS 

We have presented the results of the most detailed 
analysis of sterile neutrino dark matter production yet 
performed. In doing so, we have explored the implica- 
tions of an unexpected coincidence: the peak production 
epoch for sterile neutrino dark matter can be the same 
as the QCD epoch where the quark and gluon degrees of 
freedom disappear and give way to color singlets. The 
effects of the QCD transition on sterile neutrino relic 
densities are bracketed by simple crossovers (no phase 
transition) and first-order phase transitions for the re- 
ordering of the strongly interacting degrees of freedom 
from quarks and gluons to hadrons. 

We find that the varying alteration in the expansion 
dynamics in these different models of the QCD tran- 
sition at high temperature and low chemical potential 
can give rise to an appreciable range of closure contri- 
butions associated with a given singlet sterile neutrino 
mass, vacuum mixing with active species, and primor- 



dial lepton number (s). In particular, first-order phase 
transition scenarios for the cosmic QCD epoch produce 
a constant temperature interval with a duration of the 
order of the Hubble time associated with phase coex- 
istence. This can result in a significant enhancement of 
sterile neutrino relic densities relative to simple crossover 
transitions. 

Ultimately, we have pointed out here an unexpected 
connection between strong-interaction scale physics in 
the early universe on the one hand and low energy scale 
dark matter relic physics at the present epoch on the 
other. Whether this connection could someday be ex- 
ploited to produce insights into the cosmic QCD epoch 
would, of course, depend on the identification of a signif- 
icant component of the dark matter as being carried by 
sterile neutrinos. 

An accurate calculation of the sterile neutrino produc- 
tion entails following the evolution of the active neutrino 
scatterer number densities, sterile neutrino production, 
and entropy transfer-induced dilution for a range of pos- 
sible behaviors for the strongly interacting sector. In- 
cluding these effects, we have presented the most sophis- 
ticated calculation of the predictions of sterile neutrino 
relic densities for a range of possible QCD transition sce- 
narios and lepton number cosmologies. 
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